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ABSTRACT 

Coronal loop emission profiles are often of remarkably constant width along 
their entire lengths, contradicting expectations based on model coronal magnetic 
field strengths decreasing with height. Meanwhile Paul Bellan has produced a 
theoretical model in which an initially empty, twisted force-free loop, on being 
filled with plasma via upflow at each foot point, in the absence of significant 
gravitational effects, forms a narrow, filamentary loop of constant cross-section. 
In this paper, we focus on equilibrium states that include stratification by uniform 
gravity while retaining the effects of magnetic field twist. Comparing these with 
related force-free equilibria, it is found that injection of \ow-(3 plasma under 
coronal conditions is not likely to change the shape of a loop significantly. These 
linear equilibria apply to the interiors and boundaries of loops only, with external 
influences modeled by boundary total pressures. The effects of total pressure 
balance with surroundings and of gravitational stratification are to inhibit the 
pinching of a loop to a constant cross-section. Only if the plasma (3 were high 
enough for the plasma to reconfigure the external field and the hydrostatic scale 
height much greater than the loop size could the final state have nearly constant 
cross section. We do not expect this to occur in the corona. 

Subject headings: magnetohydro dynamics: Sun, solar magnetic fields, solar corona 



1. Introduction 

The most emissive steady solar coronal fields at MK temperatures are the closed, loop- 
shaped structures that trap coronal plasma against thermal expansion (Rosner et al. 1978, 
Aschwanden 2004). Understanding the basic physical properties of these structures is a 
crucial step towards understanding the corona. Closed plasma loops are believed to trace 
lines of magnetic force that penetrate the photosphere from below and expand to fill the 
coronal volume above an active region (Bray et al. 1991, Aschwanden 2004). They are 
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curved filamentary structures connecting regions of opposite photospheric magnetic polarity. 
Plasma loops are often observed as isolated brightened flux trajectories embedded in a less 
emissive field configuration. A broad consensus has developed that coronal loops are nearly 
isothermal along their arclengths, while temperatures and densities can vary on small scales 
between neighboring loop structures (Rosner et al. 1978, Lenz et al. 1999, Aschwanden 
et al. 1999, 2000, Aschwanden and Nitta 2000, Schmelz et al. 2001, Martens et al. 2002, 
Schmelz 2002, Winebarger et al. 2002). This can be explained by the thermal conductivity 
being much more efficient along than across the magnetic field in the nearly fully-ionized 
coronal gas (Spitzer 1962) and by the very high electrical conductivity of coronal plasma 
which "freezes" the plasma into the field lines of force (e.g. Aschwanden 2004). 

However, much remains to be learnt about coronal loops. Two physical parameters 
that continue to cause controversy are the scale height of hydrostatic stratification and the 
width of the observed loop-shaped emission patterns. Loop emission in active regions is 
frequently clearly visible at heights several times greater than the hydrostatic scale height 
associated with the temperature of emission (Aschwanden et al. 2001). Furthermore, the 
widths of these emission signatures remain remarkably constant along their entire lengths 
(Klimchuk 2000, Watko & Klimchuk 2000). These emission patterns are usually identified 
with magnetic flux tubes because of the efficient field-aligned thermal conductivity of the 
corona, and they therefore defy theoretical expectations that flux tubes must expand in a 
stratified magnetic field (Lopez Fuentes et al. 2006). In a divergence-free flux distribution 
(V • B = where B is the magnetic field vector), flux tube cross-sectional area and field 
strength must be inversely proportional. Since the macroscopic coronal magnetic field may 
be approximated by a superposition of magnetic multipoles (Altschuler & Newkirk 1969) 
and there is obviously more space high in the solar atmosphere than low down, we expect 
the field strength to decrease with increasing altitude. Pioneering direct measurements of 
coronal magnetic field strength by Lin et al. (2004) also show a marked decrease of field 
strength with increasing altitude. Yet loop emission profiles generally have near-constant 
width. While concerns have been raised about the effects of background subtraction and 
spatial resolution on the measurement of the widths of these narrow structures, in practice 
width variations and background are not correlated and loops with constant cross-section 
can be distinguished with confidence from loops with expansion factors characteristic of 
force-free models, even if these widths are close to the resolution limit (Lopez Fuentes et 
al. 2008). 

The variation of cross-sectional area along a loop is particularly important to under- 
standing plasma dynamics along loops. Klimchuk et al. (2000) showed that the tension force 
associated with twist tends to make more twisted fiux tubes more circular in shape, but 
that the near-constant cross-sections of loops cannot be explained by this twist. Petrie et 
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al. (2003), Patsourakos et al. (2004) and Gontikakis et al. (2005) recently conducted different 

studies of how steady plasma flows along loops may influence the observed fall-off of plasma 
density, and therefore emission intensity, with increasing height. Patsourakos et al. (2004) 
assumed that loop cross-sectional constant or nearly constant along loops, invoking 

the observations of Yohkoh/SXT and TRACE loop widths by Klimchuk (2000) and Watko 
& Klimchuk (2000). Traditionally HD models have assumed a constant cross-section (Or- 
lando et al. 1995a, b; Reale et al. 2000a, b, Peres 2000, Winebarger et al. 2003). Cargill & 
Priest (1980, 1982) investigated the effects of flux tube cross-sections varying hnearly with 
arclength. Motivated by interplay between characteristic scale heights of plasmas and mag- 
netic fields, Petrie (2006b) investigated steady flows of stratifled plasmas in flux tubes whose 
cross-sectional widths decrease exponentially with increasing altitude. The outcome is that 
flux tube expansion can enable the extension of plasma scale heights by steady flows, and 
such flows can extend scale heights to observed scales, but that the conditions under which 
this can occur are rather special and are unlikely to be met generally in loop structures. 

In this paper we address the issue of how loop-shaped emission patterns of constant 
cross-section are created. In particular, we investigate whether or not a mechanism proposed 
by Bellan (2003) for the creation of thin, fllamentary flux tubes of constant cross-section, 
formulated in the absence of gravity, can work in the gravitationally stratifled medium of 
the corona. Bellan presented a sequence of physical processes which result in a thin, hot 
flux tube of exactly constant cross-sectional width. The model is consistent with the thin, 
fllamentary structures of solar loop emission patterns as well as with various experiments 
performed by Bellan's plasma physics group (e.g. You et al. 2005). 

In axisymmetric geometry with cylindrical coordinates {r,9,z), where the symmetry 
axis r — coincides with the central axis of the loop flux tube, the process begins with a 
current-free loop, i.e., — 0, with held stronger at the foot points z — ±L than at the 
apex z — 0. This means that the tube is wider at the apex than at the ends initially. If a 
current / flows along this tube then a toroidal fleld component Bg ^ is introduced and the 
fleld becomes twisted. The resulting non-conservative Lorcntz (j x B) forces then accelerate 
plasma from the narrow ends of the tube towards the apex. In Bellan's (2003) picture, 
these flows thereby advect frozen-in toroidal magnetic flux from each foot point towards the 
apex so that there is excessive magnetic tension there. This tension pinches the loop until 
a new equilibrium is found. In the absence of gravity the only possible equilibrium is one- 
dimensional and independent of 2;: it has constant plasma and magnetic flux distributions 
along the entire length of the loop and so we arrive at a thin loop of constant cross-section. 

An important difference between the solar case and the case studied by Bellan (2003) 
is the presence of a signiflcant gravitational force in the solar corona. While the coronal 
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gas is not dense, the gravitational scale height of plasma at typical coronal temperatures, 
^ 1 MK, is about 60 Mm, where we assume a fully ionized hydrogen plasma. Aschwanden 
et al. (2001,2004) quote a value of about 47 Mm, assuming a 10:1 ratio of hydrogen to 
helium. Therefore, loops of active region scales, of order 100,000 km, should be modeled 
with gravitational effects taken into consideration. 

The dynamics and the final equilibrium problem described above are affected by the in- 
troduction of gravity. The acceleration of plasma from the thin parts of the loop towards the 
apex by the Lorentz forces is aided by the plasma pressure gradient force and is therefore, in 
the absence of gravity, unopposed. The addition of a gravitational force that acts downwards 
from the apex towards foot points introduces opposition to this acceleration. For even very 
slight loop twist, Lorentz forces in this case with gravity are still likely to be large enough 
for the acceleration to take place and fill the loop (Bellan 2007, private communication), but 
the manner in which the loop is filled with plasma and the subsequent equilibration will be 
modified by the gravity. Without gravity, the equilibrium equation is j x B = Vp, with j 
and p denoting the electric current density and plasma pressure, so that one can immedi- 
ately see that B ■ Vp = and an isothermal plasma must be evenly distributed along each 
magnetic flux trajectory. In the presence of a gravitational force some of which acts along 
the field direction, the plasma is stratified by gravity and cannot be uniformly distributed 
along the fields. As well as its non-uniform distribution, the possibility of a low-/? plasma to 
reshape a loop's boundary with a dynamically dominant external magnetic field must also 
be questioned. 

This equilibrium problem is the focus of this paper, which is organized as follows. The 
mathematical problem is developed in Section [2] and solutions presented in Section |3l We 
conclude with a discussion in Section HI 



2. The magnetohydrostatic problem 

Consider the static equilibrium model based on the one-fluid ideal hydromagnetic de- 
scription, denoting the magnetic field, plasma pressure and density by B, p and p, respec- 
tively. The balance of forces is described by 

-L(VxB)xB-Vp-p(7z = 0, (1) 

An 



assuming a uniform local gravity of acceleration g in the z direction, along the axis of the 
loop. Then the ideal gas law relates the gas pressure p to the gas density p 
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P = ^PT, (2) 
/i 

where is Boltzmann's constant and fi is the mean particle mass for a fully ionized 
(monatomic) hydrogen plasma. The solenoidal condition 



V • B = 0, (3) 

closes the set of equations to determine p, p, and B. To keep the physical problem simple, 
we avoid the complication of a full energy equation by applying the isothermal case T = Tq 
a constant. 

In the 2D case in cylindrical geometry with axisymmetry {d/dO = 0), the magnetic field 
satisfying equation has the form 

B = ^ f^r, z), Mllil^ , (4) 

r \ dz ' ' ' dr J ' 

where f{r,z)/r is the toroidal magnetic field component and if){r,z) is the magnetic flux 
function, whose isosurfaces contain the magnetic field trajectories. By standard theory (Low 
1975), the momentum equation ([T]) reduces to 



d'^ip 1 dip d'^ip 



r dr ^ dz"^ 



dfW , ^_dp{ip,z) 
dip 



+ 4tt 



dp{ip, z) 



dz 



dip 

+ p{'^,z)g{z) 



0, 
0. 



(5) 
(6) 



i/i=const 



Here f{r,z) = f{ip) must be a strict function of ip from the projection of the momentum 
equation on the symmetry axis. 



3. The models 

This section focuses on the static equilibria before and after a loop is filled with plasma. 
This paper does not directly treat the dynamics of the evolution to a final equilibrium, but it 
is shown that the final equilibrium can be reached from the initial equilibrium via ideal MHD 
processes. The initial equilibrium model describes an empty loop, i.e. p = p = 0, and so there 



- 6 - 



are no forces in this equilibrium. The initial equilibrium is therefore described by a force-free 
field. The final equilibrium does involve forces and requires a magnetohydrostatic (MHS) 
description. We describe linear force-free models for the initial equilibrium in Subsection 13.11 
and related linear MHS equilibria in Subsection 13.21 



3.1. Linear force-free fields 

The force-free case is the one with po(^) = constant (zero in our scenario) so that all 
Maxwell stresses are contained within the magnetic field and all electric currents are aligned 
with the magnetic field. The current associated with the magnetic field of equation (jlj) is 
described by 

j = Bo f-V'^, W) _dmdl\ ^ 
\ ^ dip dz^ dip dy J ^ 

The general condition for the nonlinear field to be force-free is that equation ([5]) is satisfied 
with dp{r, z)/dz = 0, in which case j = a(r, z)B with the nonlinear force-free parameter 
a{r,z) = df{ip)/dip - see equation ([5]). The field-aligned currents are clearly associated with 
axial magnetic fiux Bg which is responsible for twist and shear in the magnetic field in these 
2D solutions. For a force-free field to be solenoidal, a must be constant along each fiux 
trajectory: B ■ Va = 0. The simplest case is the potential case f{ip) = a = which has no 
electric currents, no azimuthal field component and therefore no twist. Of interest here is 
the linear case f{ip) = <yip, for some constant a, which is seen from equation (j5]) to yield the 
usual Helmholtz equation for linear force-free fields {'V'^ + a'^)ip = 0. For a narrow equilibrium 
structure like a loop, this linear case should be sufficient. 

Browning & Priest (1983) found three linear force-free solutions of equation ([1]), one of 
which lends itself to the modeling of an expanding coronal loop: 



ipff{r,z) = iporJi{kr) cosh.{mz), = k'^ — a'^ , (8) 

where a is constant (henceforth we will only treat constant-a cases). Three examples of the 
solution described by equation ([8]) are plotted in Figure [H for different values of m but all 
with a = 1/16. We expect that during dynamical phases, the loop will remain anchored 
to the base of the atmosphere since the photospheric plasma is orders of magnitude denser 
than the coronal plasma, and should not be affected by coronal field evolution of this type. 
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The footprint of the loop at the base of the corona z = ±100 Mm is fixed to be the circle 
r < a Mm by setting 



aJi{ka) cosh(mL) ' 

where L = 100 Mm is the loop half-length, and by adopting the surface ip = 1 as the 
loop boundary. Here we show examples with a = 4 Mm,which corresponds reasonably to 
observations by Aschwanden et al. (1999) and Watko & Klimchuk (2000) and will allow 
other physical parameters to take reasonable values as described later. In Figure [2] the top 
left picture shows ip\z=±L- the cases k = (a^ + m^)^/^ for a = 1/16 are graphed as the 
dashed lines, all almost coinciding with the solid line, and a = 1/4 (dotted lines). Each 
of these distributions is shown for three values of m: m = 1/25 Mm, m = 1/50 Mm and 
m = 1/100 Mm, corresponding to the left pictures of Figure [H However, the differences 
between these different values for m are small in the graphs because a is more dominant 
than m in the expression for k: for the length scales of interest, the twist affects the boundary 
flux distribution more than the axial variations do. For small k, ipz=±L is very close to the 
solid line, which describes the parabola i(j = r'^. This is the flux function (at z = ±L) of 
the related MHS equilibrium, which we introduce in the next subsection. The magnetic field 
vector is 



B// = _Bo [mJi{kr) sinh(mz), aJi{kr) cosh(mz), kJo{kr) cosh(mz)] / [aJi{ka) cosh(Lm)] . 

(10) 

The inverse length scales k, a and m all influence the character of the solution in 
distinctive ways. The inverse length scale m controls the shape of the loop: increasing m 
gives the loop a more of a bulging cross-section because Br is proportional to m. Figure [2l 
(top right) shows example distributions of Br\z=-L for selected cases of the force-free solution: 
k = (a^ + m^)^/^ for a = 1/16 (dashed lines) and a = 1/4 (dotted lines). These graphs 
are shown for three values of m: m = 1/25 Mm, m = 1/50 Mm and m = 1/100 Mm, 
corresponding to the left pictures of Figure [H The equivalent graphs for Br\z=+L are the 
mirror image of these in the r-axis. The correspondence between the increasing loop bulge 
and the increasing Br as m increases is clear. The a parameter, when non-zero, adds twist 
to the magnetic field, since Bg is directly proportional to a. The bottom left picture of 
Figure [2] shows example distributions of B0\z=±l with a = 1/16, the value of all force- free 
examples in Figure[TJ Shown are cases with k = (a^ + m^)^/^ for a = 1/16 (dashed lines) and 
a = 1/4 (dotted lines) for comparison. The case a = 1/16 is representative of all examples 
in Figure [1] (left pictures) since the dashed lines are close together, particularly within the 
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loop. The axial flux is proportional to the parameter k. Example distributions of Bz\z=±l 
are shown in the bottom right picture of Figure [2|, again for k = (a^ + m^)^/^ with a = 1/16 
(dashed lines) and a = 1/4 (dotted lines). The dashed lines are close to each other, and 
correspond to all three pictures in the left half of Figured! The distributions of Bz\z=±l and 
Be\z=±L are effectively independent of m while the distribution of Br\z=±L is not. This is 
because Br is related to the bulge of the loop to a significant degree and Bg and B^ are not. 

The three parameters k, a and m are not independent of each other: the toroidal flux 
is limited by the axial flux strength via the condition = k"^ — > which then fixes the 
radial field strength. A ID configuration is possible if the toroidal field strength provides 
exactly strong enough magnetic tension to balance the expansive magnetic pressure gradient 
force of the axial flux, when a = k. This is a limiting case beyond which the toroidal flux is 
too strong for an equilibrium solution to exist. Increasing k increases the expansion of the 
flux tube while increasing a draws in this expansion. The parameter k, and these ip and 
B distributions generally, will receive further attention after we have introduced the MHS 
solutions in the next subsection: for small k, the field component distributions are strikingly 
similar to the related MHS solutions that form the subject of that subsection. In practice, 
length scales of variation along loops are so large that the inverse length scale m is generally 
very small, and interesting cases have k ^ a and a > 0. At least as far as this solution class 
is concerned, electric currents are necessary for interesting physics. Bellan (2003) needed a 
toroidal field component to be present for his process to work. 

We omit here the effects of a curved loop axis, assuming that the effects of curvature are 
not important to the physics under consideration. This assumption has been applied many 
times in the study of twisted loop equilibria in cases with or without forces (Parker 1979, 
1994, Browning & Priest 1983, Zweibel & Boozer 1985, Browning & Hood 1989, Lothian 
& Hood 1989, Parker 1994, Belien et al. 1997a, b. Van der Linden & Hood 1999, Bellan 
2003). The current self- repulsion or hoop force characteristic of curved current structures 
is unlikely to be important in for the physics of a typical loop: in the scenario described 
by Gold (1964) and Parker (1994), a loop is one modestly twisted flux tube among many 
flux tubes, probably separated by current sheets in general, and a loop's self-repulsion is 
likely to be much weaker than local pressure balance between neighboring flux tubes which 
will squeeze the loop or else allow it space in which to expand. The force-free solutions 
of the previous section demonstrate the tendency of force-free fields to expand to fill an 
atmosphere, unless reined in by the tension force of strong azimuthal flux. Our objective 
here is to determine whether or not plasma forces increase the likelihood of a slightly twisted 
flux tube assuming a constant cross-sectional shape in a gravitationally stratified medium, 
as Bellan (2003) found in the case without gravity. 
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Fig. 1. — Linear force-free solutions (left pictures) and magnetohydrostatic (MHS) solutions 
including forces (right pictures) corresponding to stratification length scales (1/m for force- 
free solutions, A for MHS solutions) of 25 Mm (top), 50 Mm (middle) and 100 Mm (bottom). 
All solutions plotted here have a = 1/16 Mm~^. 
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Fig. 2. — A direct comparison between the flux function (top left) and magnetic field com- 
ponent distributions (the three other pictures) at the boundary z = —L of the MHS solution 
(solid lines) with various cases of the force-free solution: k = (a^ + m^)^/^ for a = 1/16 
(dashed lines) and a = 1/4 (dotted lines). The force- free and MHS solutions are plotted for 
the cases m = 1/A = 1/25 Mm~^, 1/50 Mm~^ and 1/100 Mm~^. In the top right picture, 
these cases are distinct because of the proportionality of Br to m. In the other pictures, the 
force-free cases coincide to a good approximation, while the MHS cases are exactly equal . 
For comparison, all force-free solutions plotted in Figure [1] have a = 1/16 Mm~^. 
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Fig. 3. — Shown by the dotted hne is the magnetic pressure (and therefore total pressure) 
on the loop boundary = 1 aX the apex 2; = as a function of hydrostatic scale height A for 
our special class of MHS equilibria. This apex boundary pressure increases as A increases 
and the loop solution becomes narrower. As A — > 00 this pressure tends to a constant, 
represented by the dotted line. This is the boundary pressure of Bellan's (2003) solution, 
the narrowest possible solution. The total boundary pressures for the equivalent force-free 
solutions with m = 1/A are smaller. Therefore, the much enhanced boundary pressure for 
large A makes it impossible for a loop to evolve from an expanded force-free configuration 
to a non-force-free configuration with significantly reduced expansion. 
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3.2. Inclusion of forces associated with isothermal plasma 

We now turn our attention to non-force-free MHS solutions. Equation ([6]) is solved by 

p{ij,z)=pM- [ p{^,z')g{z)dz'. (11) 

In general, the Grad-Shafranov equation ([5]) can assume as many different forms as there are 
distinct regions of physical space, each of which may be occupied by a different branch of the 
functions p{tp, z) and p^ip, z). The distinct regions are separated by free boundaries, the loci 
of which are unknown and generally can only be calculated simultaneously with the solution 
if). The free boundaries are governed by continuity or jump conditions imposed by force 
balance at these boundaries. Free boundaries also separate regions where the functional 
forms of p{iIj,z) and p{ilj,z) differ, e.g., a region of potential or force-free magnetic field 
separated from a region of full magnetohydrostatic force balance. The nonlinear equations 
involved require a full numerical treatment along the lines of Petrie et al. (2007). We defer 
such a study and choose to stay close to the linear treatments of past authors, including 
Bellan (2003), to enable a direct comparison. 

For simplicity we set P{ip) = Aip, so that equation ([5]) is linear. In the isothermal case, 
the hydrostatic scale height A = p/ (pg) is constant. Now the condition p/p = constant gives 
PiW = in equation (ITT]) . 

We seek only separable forms for p{tp, z): 



p{ij,z)=poWs{z). (12) 
The Grad-Shafranov equation ([5]) becomes 

Qj.2 ^ Qj, Q^2 

Equation (fT3|) generally needs to be solved numerically (Belien 1997a, b). However, if we 
prescribe g{z) using an elementary analytical function: 

^(z) = ^otanh(^|) , (14) 



where go is the solar surface gravitational acceleration 274 m/s^, then g{z) models the effect 
of gravity acting in opposite directions along each leg of the loop, in the negative direction 
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in the left leg [z < 0) and in the positive direction in the right {z > 0) leg. This form of 
g{z) tends to the uniform surface gravitational acceleration at the loop foot points z = ±1, 
where the field is assumed to be nearly vertical, and changes sign at the apex of the loop 
z = where the loop is locally horizontal. This g{z) therefore captures the global effect 
of gravitational stratification on a twisted loop while ignoring the component of gravity 
perpendicular to the curved loop axis. This is reasonable for the loops of large aspect ratio 
(length/ width) that concern us here, since the effects of axis curvature are not large for such 
loops (Belien et al. 1997a, b). 

Now equations flll|14p together determine the dependence on z of p. With g{z) given 
by equation (|T^ . s{z) must take the form 

siz) = cosh (I) , (15) 

which, in the near-vertical portion of the loop near the foot points models the exponential 
fall-off of pressure and density characteristic of hydrostatic gravitational stratification. The 
pressure and density are both minimized at the apex 2; = 0. A further advantage of this form 
for g{z) is that equation ( fT3l) becomes separable such that exact analytical solutions can be 
obtained, and the resulting MHS solution can be easily compared to the linear force-free 
field solution of equation ([8]) with m = 1/A. Defining ip = tpi a.s the flux surface on which 
p vanishes and scaling ip so that = 1 on this surface, if we assume also a linear form for 



PoW=Po{l-^), (16) 
then we have reduced the problem to a linear equation: 

ITT - -IT + ITT + ^ -^'^^or =0, (17) 
ar^ r or oz'^ 

We seek a particular integral of the Grad-Shafranov equation (ITTI) that balances the 
source term AnP^r"^. If we fix -j/^ = 1 on the flux surface passing through r = a at 2; = ±L 
where L = 100 Mm is the loop half-length, then 



0^2 



cosh(L/A) 



where L = 100 Mm is the loop half-length. This particular integral is therefore 
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if^mhs = r'^ cosh I cosh f ^ ) , (19) 



and the magnetic vector is 



— sinh ( ^ ) , ar cosh ( — ) , 2 cosh ( 4 , , - , . 

A VA/ VA/ VA^J / V A 



a' cosh [\\ . (20) 



We plot example solutions in Figured] (right pictures). The effect of gravitational stratifi- 
cation is clear: the plasma is concentrated close to the foot points of the loop and the loop 
is not of uniform cross-sectional width. We have chosen here cases with hydrostatic scale 
height A = 25, 50 and 100 Mm corresponding to temperatures of about 0.5, 1 and 2 MK, 
respectively. In the limiting case with A ^ oo, corresponding to a very hot plasma (or 
a scaled-down loop) we recover the axially uniform case found by Bellan (2003) in which 
gravity does not play a significant role. 

The plasma (5 on the loop axis is given by 

+ (21) 



Bl y^^o 2 VA' ; cosh(L/A) 

which at the boundaries z = ±L is equal to (l/A^ + a')/2 ^ 0.044, 0.034 and 0.032 in the 
three cases A = L, L/2 and L/4, respectively in Figure [1] with a = 4 Mm. At the apex 
2; = 0, /? is a factor of cosh(L/A) smaller, and falls to about 0.021, 0.0092 and 0.0016 for 
the three cases A = L, L/2 and -L/4, respectively. Obviously these /5 values decrease with 
distance from the loop axis r = 0. 

The angle between a field line and the axial direction at any point, the pitch angle, is 
where tan(0) = Bq/Bz- A field line makes one complete revolution around the loop axis in 
a distance 

, / Bz 2'Kr , , 

\{r,Q) = 2Tir-^ = -. 22 

^ ' Bo ta.Yi(h ^ ' 



Thus 



A 



2TTrkjQ{kr)/aJi{kr), force-free case, 
Air/a, MHS case. 
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In all of our examples, a = 1/16 Mm~^, and so in our MHS solutions A = QAn Mm and the 
lines turns through nearly 27c radians along the entire length of each loop. This is within 
the Kruskal-Shafranov limit of 2tt turns for a simple toroidal configuration (Kruskal 1954, 
Shafranov 1957), often considered a representative condition for stability against kink modes 
(e.g. Priest 1982, Goedbloed & Poedts 2004). In the force-free models considered here, with 
A; ~ a = 1/16 Mm~^, the result is similar. The stability properties of these solutions can 
be determined in the form of a full resistive MHD spectrum by solving the linearized MHD 
equations using the hyperbolic stability solver PHOENIX (Blokland et al. 2007). This will 
be the subject of a later study. 

This set of non-force-free MHS solutions does not include the force-free solution of 
equation ([H]) as a special case: setting Pq = in equation f[T^ . implying that 1/A^ + = 0, 
just gives the trivial solution ip = 0. However, the two solutions can still be related to each 
other in a meaningful way. Compared to this non-force-free solution, the force-free solution's 
photospheric axial flux distribution is more concentrated towards the central axis. For small 
k the force-free solution is very close to the non-force-free solution but, since k = forces 
m = a = 0, k>0 for any non-trivial solution. The resemblance between the force-free 
and non-force-free solutions plotted in Figure [T] is due to the small size of k. Since axial 
length scales in the solar corona are so large, measured in tens of Mm, this means that 
is small and so k and a are of similar size in the examples plotted in Figure [TJ In the 
top left picture Figure [2] the magnetic flux functions ipff and ipmhs, given by Equations ([9]) 
and flT^ .are graphed for comparison of the non- force- free and force-free solutions. While 
these graphs show photospheric distributions, multiplying by cosh(mz) or cosh(z/A) gives 
the force-free or MHS distribution for any —L < z < L. For A; < 1, Ji{kr) is nearly linear in 
the domain of interest, and so ipff and ipmhs are nearly equal. Since Ji has positive second 
derivative in this domain of interest, ipff is slightly larger than ipmhs everywhere between 
r = and r = a except on the axis r = where both are zero, and at the outer boundary 
where both solutions are 1. This implies that evolving from ipff to ipmhs would involve a 
small outward movement of each flux surface ip = constant wherever r < a, as expected 
when gas is pumped into a flux tube adding to the outward pressure force. On the other 
hand, for r > a, ipff is smaller than ipmhs, as indicated by Figure [2], meaning that all flux 
surfaces outside the cylinder r = a actually move inward during evolution from tpff to ipmhs- 
This is because there is less plasma out here than close to the axis, so that outward pressure 
forces are not greatly enhanced here, and the outward displacement of axial magnetic flux 
by the plasma causes the axial fleld component to be stronger than in the force-free case, 
without enhancing the radial component, so that the flux tube has less bulge. This effect is 
greater for greater m and smaller A, in which case more of the tube lies outside the cylinder 
r = a and is therefore subject to this inward movement. This difference between empty and 
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full flux tubes is therefore reminiscent of the picture presented by Bellan (2003) of plasma 
injection causing a flux tube to become thinner, although the change is much smaller than 
in the case focused on by Bellan. 

All of the MHS examples of Figure [1] (right pictures) have Bz\z=±l = Bq/8 G and 
Be = Bor /256 G. These are the solid curves in the bottom pictures of Figure [2l Within 
the loop the boundary distributions of the axial and azimuthal fleld components Bz\z=±l 
and B0\z=±L are very similar for the force-free models represented by the dashed curve in 
each picture, and the MHS models represented by the solid curve. The force-free axial fleld 
distribution decreases very slowly with radius, is larger than the MHS value at the center of 
the loop and crosses the MHS value aX r ^ 0.7a Mm. Meanwhile the azimuthal component 
is slightly greater than the MHS component for < r < a, matching it at r = 0, a. 

Unlike Bg\z=±L and Bz\z=±l, Br\z=±L is strongly affected by changes in m or A. Figure [2] 
(top right picture) shows the radial fleld distributions at z = —L for the three values of 1/m 
and A represented in Figure [TJ 25 Mm, 50Mm and 100 Mm. Br is proportional to this 
parameter and so the greater the bulge of the loop, the greater the radial fleld component. 
Within the loop the MHS radial fleld of a given value of A = Aq is slightly smaller than the 
force-free radial fleld with m = I/Aq except at r = 0, a Mm where they match. This agrees 
with the MHS loops being slightly less bulged than the force-free loops. 

In the force-free solution class, m depends on the twist parameter a and the radial flux 
distribution associated with k (Figure [2]), whereas in the non-force-free solution class A is 
determined entirely by the temperature of the injected plasma. Since these two parameters 
are independent, any of the force-free solutions plotted in Figure [1] could evolve to any of 
the non-force-free solutions, provided that certain invariants of the ideal MHD equations 
are conserved and certain boundary conditions are met. In the corona, conditions nearly 
correspond to those of ideal MHD, meaning that evolution from one static state to another 
must preserve the total flux through the tube, as well as the helicity trapped in the magnetic 
fleld to a very good approximation. The latter condition is met by imposing the same 
constant value of a on beginning and end states. The former is guaranteed by the fact that 
we have normalized each solution so that the total flux is equal in all solutions in Figure [TJ 

Not only must the total magnetic flux and helicity in the loop be conserved, but the 
boundary distribution of vertical flux Bz\z=±l and the distribution of twist Bg\z=±L must be 
approximately conserved also. These conditions are physically imposed on the coronal fleld 
by the inertia of the heavy plasma in layers beneath, which forbids the axial flux and twist 
from being signiflcantly redistributed. As we have discussed. Figure [2] (bottom pictures) 
shows that these distributions are approximately equal in the vicinity of the loops for all 
examples, force-free and MHS, shown in Figure [T|, and that there is no difference at all 
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between the distributions for the force-free examples shown in the left pictures of Figure [H 
nor are there any differences between the these distributions for the MHS examples of the 
right pictures. In the domain ^Z' < 1, where both classes of model are valid, the force- free 
and MHS distributions of Bz\z=±l and Bg\z=±L are approximately the same for the modest 
values of ~ a expected in the corona. 

Therefore, under the restrictions described so far, each solution in Figure [1] could in 
principle be ideally evolved into any other. On the other hand, an important ingredient 
missing from these models is a suitable external magnetic field in the MHS models. Those 
presented here have unphysical plasma pressures that distort the external field - an effect 
that grows with distance from the loop. Consideration of our MHS magnetic fields outside 
the boundary flux contour ip = 1 is not strictly meaningful because for ip > 1, in view of the 
linear form of Equation flTB]) . the plasma pressure p becomes negative, a deficiency shared 
with the few existing MHS loop flux tube models (Belien et al. 1997a, b, Bellan 2003) whose 
expressions for the plasma pressure take a linear form. Likewise, the force-free external fields 
are not useful as models of coronal phenomena outside the loop. In Figure [2], the unphysical 
a = 1/4 curves for < r < 20 Mm give a good idea of what our a = 1/16 curves look 
like in < r < 80 Mm. In MHS and force-free cases we cannot expect the distant coronal 
field to wind around a single central axis. A model combining a loop in MHS or twisted 
force-free equilibrium with a reasonable force-free external field will necessitate a numerical 
treatment to cope with the nonlinearities involved, along the lines of the prominence models 
of Petrie et al. (2007). In the meantime, we can interpret the analytical solutions presented 
here, taking into account the likely influence of the external coronal field. 

At the boundary ip = 1 the loop is constrained by the external field, near-potential 
fields or other twisted loops not included in the models but whose effects on the loop can 
still be anticipated. In reality, while the field inside a loop is relaxing to equilibrium with 
its newly injected plasma, the force-free field around the loop will react to the changing 
shape of the loop, moving to fill space vacated by the loop if it narrows, or else being pushed 
outward by the loop if it widens. Although the external fields in these linear solutions are not 
appropriate to the physics of a real external coronal field and plasma, the essential physics 
can be inferred from the total pressure imposed on the loop by the external field at the 
boundary. In particular, at the apex z = the total pressure on our MHS loop model takes 
a simple form dependent only on A: 



lp.|2 _ p2 4 + aV cosh(L/A) 



This quantity is graphed as a function of A with a = 1/16 Mm, a = 4 Mm and L = 100 Mm 
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in Figure [31 Clearly this pressure increases significantly as the temperature and the scale 
height A increases and as the model loop becomes narrower. The extreme case A ^ oo, 
Bellan's (2003) model, has a much higher boundary apex pressure, represented by the dotted 
line, than any of the models in Figures [T] and [H On the other hand. Figure [2] shows that 
the equivalent force- free fields with m = 1/A have smaller boundary pressure than their 
MHS counterparts. This means if that a force-free solution in a given row of Figure [1] 
with a given value of m is injected with plasma with scale height A = 1/m, the final MHS 
equilibrium cannot be as narrow as the one on the same row of Figured] because its boundary 
pressure is too strong. If the loop were to become so narrow, there would be a total pressure 
imbalance at the boundary, and so the boundary would move outward until pressure balance 
were restored. In the same way, if plasma of arbitrary scale height A > 1/m is injected, the 
correspondingly narrower equilibrium could not be attained because the total loop boundary 
pressure would be too large for the external field total pressure. Thus it is clear that an 
expanded force-free loop in Figur^H on being injected with low-/? plasma, must attain final 
equilibrium with higher expansion than either its MHS counterpart in Figure [1] or the MHS 
solution corresponding to the scale height A of the injected plasma, neither of which solutions 
is of uniform cross section at EUV temperatures. If the MHS solution corresponding to A 
is much narrower than the initial force-free solution then the actual final MHS state must 
have much smaller total boundary pressure than this solution and must therefore be much 
more expanded. Therefore we can only conclude that plasma injection cannot explain the 
uniform cross-sectional profiles that we see in images of loops. 



4. Discussion 



We have revisited the long-standing mystery of why loop-shaped coronal emission signa- 
tures have constant cross-sectional width, while models lead us to expect coronal flux tubes 
to have width increase with height. This work has concentrated on the effects of a uniform 
gravitational field on the magnetic field and plasma distribution of a moderately twisted 
coronal loop. Bellan (2003) argued that an initially empty twisted loop, on being injected 
with plasma in the absence of gravity, must evolve to a thin equilibrium state of constant 
cross-section. We examined equilibrium state including effects of gravitational stratifica- 
tion, motivated by the fact that measurements of widths of EUV loops at about 1 MK are 
consistent with constant cross section even though their hydrostatic stratification scales are 
comparable to or smaller than the sizes of the loops themselves. Introducing gravity into 
the modeling results in equilibria, with MHS forces, of non-constant cross section in gen- 
eral, with near-constant cross section only when the loop is hot enough for the stratification 
length scale to be significantly larger than the size of the loop. Taking into account the 
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influence of the external coronal field, the injection of plasma into an empty, expanded loop 
could not result in the loop becoming significantly narrower, certainly not narrow enough to 
have nearly constant cross-section. Therefore, while we have not solved the mystery of the 
constant-width loop, we have decisively eliminated plasma forces as a possible cause for this 
phenomenon for large (> 100,000 km) at typical coronal (~ 1 MK) temperatures. 

The MHS solutions applied here arc restricted to a special linear class of solutions, which 
correspond to a specific radial distribution of injected plasma. Generally plasma injection 
into one of these linear force-free loops would result in a highly nonlinear MHS state not 
describable in closed analytical form. Also we neglected the influence of surrounding flclds, 
following many previous authors in taking the loop boundary as the boundary of the problem 
itself. The value of the special linear MHS solutions applied here hes in the clarity with which 
they bring out the essential physical consequences of the plasma flows into the force- free loops 
and of the plasma forces on the final equilibrium. 

The theoretical basis for expecting coronal flux tubes to expand is clear. Near-constant 
flux tube cross sections are inconsistent with solenoidal ( V • B = 0) field strengths decreasing 
with height. Field strength and flux tube width are inversely proportional for any solenoidal 
field, and so a field whose strength decreases with increasing altitude should have a corre- 
sponding increase in flux tube width. A dipole fleld strength falls off with increasing radius 
as while higher-order potential flelds including active-region loops should fall off more 
steeply still. Thus even taking moderate macroscopic nonlinear currents into account, loop 
widths should generally increase with height as quickly as at least. 

This leads us to question whether the loop-shaped emission patterns of constant width 
really represent flux tubes. Perhaps these patterns delineate not flux tubes but loci of mag- 
netic fleld diffusions that, locally to very narrow regions, heat successive narrow loop strands 
and accelerate them from these regions during the fleld reconnections (Petrie 2006). Dis- 
sipations of spontaneous current sheets are expected to occur along flux tube boundaries 
(Parker 1994) and may well trace out magnetic field trajectories. If diffusion region sizes, 
accelerations due to field diffusions, and plasma cooling times are uncorrelated with height 
in the corona, then diffusion processes should produce emission patterns of approximately 
constant width. It is hoped that coronal magnetic fleld measurements (Lin et al. 2004; 
Tomczyk et al. 2004) will attain sufficiently high spatial resolution for us to check directly 
whether field strengths within loops whose emission profiles have constant cross-section de- 
crease significantly with height. 

There are many difficulties in determining whether emission patterns really do represent 
flux tubes and, if so, of measuring the typical cross-sectional area variation with length of 
these flux tubes. Furthermore a loop is highly unlikely to have a simple cross-sectional shape 
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in the solar atmosphere or to maintain a single shape along its entire length (Parker 1994, 
Gudiksen & Nordlund 2005, Lopez Fuentes et al. 2006) and this may greatly influence its 
appearance from certain angles. The identification and characterization of individual loops 
will hopefully be aided by NASA's STEREO mission. 

I thank the referee for helpful and constructive comments, and in particular for stressing 
the importance of the effects of external fields. 1 also thank Paul Bellan and Eve Stenson for 
illuminating discussions. Some of this work was conducted while the author was a participant 
in the National Aeronautics and Space Administration (NASA) Postdoctoral Program at 
Goddard Space Fhght Center, and was based at National Solar Observatory, Tucson. 
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